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Studies of Feshbach resonance phenomena in fermionic alkali gases have drawn heavily on the 
intuition afforded by a Fermi-Bose theory which presents the Feshbach molecule as a featureless 
Bose particle. While this model may provide a suitable platform to explore the ^Li system, we argue 
that its application to '"^K, where the hyperfine structure is inverted, is inappropriate. Introducing a 
three-state Fermi model, where a spin state is shared by the open and closed channel states, we show 
that effects of "Pauli blocking" appear in the internal structure of the condensate wave function. 
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Fermionic alkali atomic gases present a unique environ- 
ment in which to control and explore the crossover be- 
tween BCS and Bose-Einstein condensation (BEC) 
Already the creation of a molecular BEC phase from a de- 
generate Fermi gas of atoms has been reported by several 
experimental groups Q, while studies of fermionic pair 
condensation in the crossover regime are under way 0. 
The facility to control the strength of the atomic pair 
interaction in the Fermi system relies on a magnetically- 
tuned Feshbach resonance (FR) phenomena involving the 
multiple scattering of atoms from open channel states 
into a molecular bound state formed from neighboring 
closed channel states. Current theories of the FR treat 
the molecular bound state as a featureless bosonic par- 
ticle, and characterize the total system by a Fermi-Bose 
theory familiar from studies of polariton condensa- 
tion ly| as well as models of bipolaronic superconductiv- 
ity |3- While the Feshbach molecule (FM) involves spin 
states different from the scattering states, the molecular 
boson can be regarded as distinct. However, if a spin 
state is shared, the validity of the Fermi-Bose theory as 
a microscopic model of the FR is called into question Q . 
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FIG. 1: Atomic states involved in the FR for *°K, where the 
Fermi gas is initially prepared in the two lowest eigenstates. 
The coupling between states allowed by the selection rules is 
represented by a dotted line. The FM is formed from ||, — |) 



Nowhere is this scenario illustrated more clearly than 
"'"K. To understand why, let us consider the Hamiltonian 
of a single fermionic alkali atom of integer nuclear spin / 
and electron spin s = ^: 



As-I + B- (2 /^c s - I) 



(1) 



and the lowest eigenstate 
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Here A denotes the strength of the hyperfine interac- 
tion and B the magnetic field, while fi^ and /in denote 
the electron and nuclear magnetic moments respectively. 
The Hamiltonian preserves only the quantum number 
niF = iTis + mi, but the eigenstates can be labelled by 
their total atomic spin at zero magnetic field \F,mp) 
since the energy varies smoothly with field. In the ^Li 
system (/ = 1), the hyperfine interaction is positive, and 
the lowest energy states form a doublet with total spin 
F = i. By contrast, in the "^^K system (/ = 4), the 
hyperfine interaction is negative and the hyperfine struc- 
ture is inverted such that the lowest eigenstate is the one 
of highest weight, viz. F^^x = —mp = | 0- Now, 
if we ignore inelastic collisions or interactions that in- 
volve spin- flips, the interatomic interaction is specified 
by a two-body potential that depends only on the elec- 
tron spin: 

F(ri - rs) = yc(ri - r^) + K(ri - Y2) Si • S2 . (2) 

Therefore, it preserves the total spin projection of the 
two-body system Mp = mp^i + mp^2 and any scatter- 
ing process between atomic states that conserves Mp is 
allowed. If one considers only low-energy, s-wave scat- 
tering - the regime relevant to experiment - interactions 
involving identical fermions are forbidden and the sub- 
space of interacting atomic states is further restricted. 
Specifically, in the ^Li system, the interaction provides 
a mechanism to affect a FR through the coupling of the 
lowest two F = ^ (open channel) states to the higher 
energy bound state formed from all pairs of hyperfine 
states, involving the F = ^ (closed channel) states, that 
satisfy the condition Mp = \ — \ = ^- Crucially, the 
constraint on Mp is even more restrictive in ^"K allow- 
ing the two states If, — |) and ||, — |) that constitute the 
open channel to couple to only one closed channel state 



2 



(Fig.^. Thus, the FM involves only a hybridiza- 
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tion of states |f,— |) and ||,— |) which competes with 
the pairing of the scattering states ||,— |) and ||,— |). 
The aim of the present paper is to explore the integrity 
of FR phenomena in the i/iree-state Fermi system and 
assess the extent to which the nature of the bound state 
impinges on the mean-field characteristics of the system. 

Although, in the three-state basis, the majority of ma- 
trix elements of the two-body pair interaction (j^J remain 
non-zero, the low-energy properties of the system may 
be characterized by just a subset of elements. Labelling 
the spin states ||,— |), and ||,— f) by indices 

i = 1, 2 and 3 respectively, the FM is created by the 
direct density interaction U between species 2 and 3. At 
the same time, the exchange contribution g, which allows 
a transfer of particles between states 1 and 2, induces an 
effective pair interaction in the open channel. As such, 
any direct density interaction between species 1 and 3 
(repulsive in the physical system) can be subsumed into 
this contribution. Therefore, at its simplest level, the 
FR of the three-state Fermi system can be modelled by 
the Hamiltonian, 



H - ^ fJ'iNi ^ ^(eici - /"O al^CLki (3) 

i—l kz 

+ ^ C^q ak2ak'3ak'-q3ak+q2 
k,k',q 

X] [^1 akiak'3"k'-q3«k+q2 + h- 



k,k',q 

where the fermion operator a^i indexes species i, ^ 
Sk'^L'^kj and, defining Ei as the corresponding eigen- 
value of the atomic interaction (J^l, Cki = fr'li?' /2m -I- Ei. 
Since the system is not in chemical equilibrium, and 
the Hamiltonian separately conserves the particle num- 
ber and A^i 4- the free energy is characterized 
by two chemical potentials /is and = ^2 = Mi2- 
Anticipating that the coupled system is prepared with 
a roughly equal population of open channel states, we 
will use the chemical potentials to impose the condition 
Ni+ N2 = Ny, = N/2. Without loss of generality, one 
can absorb Ei and E^ into a redefinition of the respec- 
tive chemical potentials, while the detuning E2 = > 
can be used to adjust the relative energy level separation 
of state 2. Finally, for simplicity, we consider the case 
where gg = jUq. 

In the following, we will present the results of a numer- 
ical mean-field analysis of the Hamiltonian across the 
FR . However, before doing so, it will be instructive to 
anticipate some qualitative aspects of the phenomenology 
that emerge from the numerics. In contrast to the Fermi- 
Bose model, the FR Hamiltonian is complicated by 
the three-fermion character of the system, but the bare 
interaction of particles in the open channel can still be 
enhanced by the formation of a two-body resonance out 




FIG. 2: (Color online) Phase diagram of the FR Hamilto- 
nian @. The solid line shows the boundary separating the 
BEC and BCS-like phases in the dilute system as inferred 
from the variational analysis @ with uo = UoN{Eo) = 3.76. 
The points marked on the curve are obtained from the numeri- 
cal mean-field analysis in the limit of low density and, in order 
of increasing 7, they correspond to the ratios A^^i /N3 ~ 0, 30, 
73, and 90%, respectively. The intersection of the curve with 
the iz-axis translates into the binding energy of the molecular 
state associated with the bare potential Uq- The density dis- 
tributions displayed in Fig.|3|are drawn from the range shown 
by crosses at 7 = 0.1. Inset: The dependence of the scattering 
length a on the detuning v, as inferred from the numerics, can 
be well-approximated by the relation {kFa)Eo/{uc — 1^) ~ 35. 



of the three-state basis. In practice, this is achieved by 
affecting an 'optimal' rearrangement of the basis states 
wherein, by exploiting the exchange interaction, states 1 
and 2 hybridize into the orthogonal combination, 

^ici' = cos ^kflki + sin 0kak2 ' 
bl^2' = ^ sin ^kflki + cos (/ika]^^ , 

such that the condensation energy associated with the 
pairing of states 1' and 3 is maximized. In this case, im- 
posing the particle number constraint, one can propose 
the variational Ansatz for the ground state wave func- 
tion. 



m = 



n[ 



COS t)u + sin 6'k flks&^ki 
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(4) 



the integrity of which is supported by the numerical anal- 
ysis below. Here, 9k encodes the overall strength of the 
condensate while (j)k defines its distribution between the 
two pairing channels: since the open channel state 3 par- 
ticipates in both condensate fractions, (a^ai) and (0302), 
there is an inherent frustration due to Pauli exclusion not 
present in the Fermi-Bose system. Since the exchange in- 
teraction contributes indirectly to pair formation, the hy- 
bridization (as reflected through c/jk) will, itself, depend 
on the strength of the condensate. To maintain contact 
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with the physical system, we will hereafter limit our con- 
siderations to situations in which the Fermi energy of 
the unperturbed system, ep = hPkp/'^m, lies far enough 
below v that the auxiliary state 2' remains unpopulated 
in the ground state. In this case, the particle number 
constraint translates to the condition /ii2 = /i3 = /i. 

Considerable insight can be gained from analytical so- 
lutions of the variational mean-field equations in the 
dilute (BEC) and dense (BCS) Hmits (cf. Ref. Q). 
When characterized by a local contact potential U{r) = 
— UoL^6{r), such an analysis reveals a phase diagram 
characterized by three dimensionless parameters, uq = 
UqN^Eq), 7 and ly/Eg where Eq = h^k'^/2m represents 
the UV cut-off set by the range of the interaction l/feg, 
and iV(e) denotes the density of states. At low densities 
— * 0, the system develops a molecular bound state 
and enters a BEC phase when v < Vc where, defining 
f{z) = 1 — yiarctan(l/y^). 



Once again, with ep <^ v, a is determined by ^ while 



/( 



2Eo' uoii^UQ + l) 



(5) 



(see Fig. EJ. In particular, one may note that the ex- 
change contribution 7 enhances the bare interaction mq 
expanding the domain of the BEC phase while, in the ab- 
sence of a direct interaction, mq = 0, the exchange can, 
by itself, induce pairing in the open channel. 

Defining the anomalous (normal) density, Kk,ji = 
($|a_kjakil*) (Pk,jj = ($|a[-akj|^')), when deep within 
the BEC phase v <^ Vc, linearization of the variational 
equations shows that the total condensate wave function 
involves the coherent superposition of components 



Kk.i3 =- sm26'kCOS0k ^ TT? T 

2 2(eki - ^) 

1 . . , (a"i +7"')aAi3 
Kk.23 =-sm26'kSm0k ~ — — ^— 

2 2(eki - n) + iy 



(6) 



where, to leading order, the condensate order param- 
eter A13 = 7J7o X]k° '^k.ia (and the partner A23 = 
Uo X^k" '*k,23) remain unspecified. Here, for \fi\ <C i', the 
chemical potential, ^ = — 1^| (which asymptotes to half 
the molecular bound state energy), is determined by the 
self-consistency condition a^^ = 7Uo/(|/i|/i?o) with the 
coefficient a = A23/A13 determined by the relation. 



1 

Uq 



(7) 



Conversely, deep within the BCS-likc phase, for |k| < kp, 
0k — {a~^ +j~^)aAi3 cot6'k <C 1 and the condensate 
wave function acquires the familiar form 

1 . 1 aAi3 

Kk,13 - 7; Sm^^k ^ T^Jl X2 I I A |2U72 ' 

2 2((ek-A*) +l"Ai3|^)V^ 

with fi ~ ep, while Kk,23 — ('/'k/2) sin 20k- For |k| 3> kp, 
the solution converges to the low-density asymptotic 




Ai3 = ^- exp 



From the variational analysis, two striking features 
emerge: firstly, in both BEC and BCS-like phases, the 
condensate wave function is characterized by two length 
scales. Deep within the BEC regime, the FM has a size 
^063 = [£'o/(!^/2 + ImI)]^^^ while that of the molecule 
formed from open channel states, fcoCi3 = {Eo/\fJ-\y^^, di- 
verges at the crossover. In the BCS-like phase, the FM is 
increased in size A:o^23 = [£'o/('^/2 — ep)]^^'^, while the 
range of the Cooper pair of open channel states is set by 
the coherence length ^13 = 'i;i7/|Q;Ai3|. Secondly, in the 
BCS-like phase, Pauli exclusion has the effect of substan- 
tially depleting the normal density pk,22 = sin^ 6^ sin^ (pk 
and, with it, the condensate fraction Kk,23 iu the range 
|k| < kp. Both features are clearly visible in the numer- 
ically inferred density distributions below (Fig.l^J. 

With this background, let us turn to the results of the 
numerical mean-field analysis. Specifically, the ground 
state wave function |$) of the three-state Fermi system is 
determined by minimizing the free energy ($|if — ^iV|<f>) 
using a generalized Bogoliubov-Valatin transformation 
flkj = Z]^=i(wk,y/3ki + K^ijf^-kj)- Here, we take the 
most general Ansatz for the ground state wave func- 
tion compatible with the formation of a condensate, i.e. 
all elements of the matrix coefficients Mk and Wk are al- 
lowed to acquire non-zero expectation values. For con- 
venience, we choose a model potential J7q that possesses 
only one bound state (although, in the quasi-equilibrium 
system, the presence of multiple bound states will not 
change the conclusions qualitatively). We set U{y) = 
— J7o exp [— (fcor)^/2], where the range of the pair inter- 
action is chosen to be much smaller than the average 
particle separation, viz. N/{koL)^ ^ 1. 

The numerical procedure involves the minimization 
of the free energy with respect to the normal and 



anomalous densities, pk 
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and Kk,ji = 



'^k jm'^k,im where, in the s-wave approximation, the 
Bogoliubov matrix coefRcients Uk and Vk, as well as the 
densities, depend only on fc = |k|. We obtain non- 
zero values of the off-diagonal component of the density 
matrix pk,i2 which is consistent with the hybrid char- 
acter of the ground state, while the observed relations 
(^I^L'foki'l^) = /Ofc,33 and ($|6j,2,6k2'|*) = confirm 
the validity of the particular variational Ansatz |0J. For 
completeness, we note that, once ep becomes comparable 
with the detuning, the population of level 2' requires an 
adjustment of the chemical potentials 7^ ^3 to com- 
ply with the particle number constraint. In this range, 
the ground state is eventually no longer encompassed by 
the reduced variational Ansatz (0J. 

The nature of the ground state can be characterized by 
monitoring the normal density Pfc,33 and the components 
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FIG. 3: (Color online) Density distribution of (a) pu.ss and 
(b) ACfe^ia and Kfc_23 for the range of scattering lengths (kFo)^^ 
shown by the crosses in Fig.|5| At A: = 0, we have acis > K23- 
The inset in (a) shows the ratio of particles Ni/Nz in the 
ground state as a function of the scattering length (kpa)"^. 
Note that the relative weight of the 1 state on the 'BCS 
side' of the resonance increases dramatically from 30% at the 
crossover to almost 100% as (kpa)'^ — > — oo. Figure (b) inset 
shows the condensate fractions A13 and A23 as a function of 
the scattering length {kpa)~^ . 

of the condensate wave function Kfc,23, ^^^13. As in single- 
channel theories involving only two species of fermions, 
the momentum distribution interpolates smoothly from 
a BCS-like distribution at {kpa)~^ ^ — 1 to a molecu- 
lar condensate wave function in the BEC regime when 
(kpa)^^ » 1, where (kpa)^^ denotes the (inverse) scat- 
tering length (Fig. . As expected from the variational 
analysis, a key feature of the condensate wave function 
is the presence of a robust tail at high momenta which 
persists into the BCS-like phase (Fig.lSJa). (Note that, to 
infer the total occupation density, the distribution must 
be weighted by the density of states ^ k^ leading to a 
significant amplification of the tail.) The existence of 
two correlation lengths and the effects of exclusion arc 
also emphasized in the variation of the condensate wave 
function. Pauli exclusion thus enhances the population 
of a quasi-molecular component to the condensate on the 
BCS side of the transition, beyond that expected in the 
single channel or Bose-Fermi models. This is amenable to 
rough estimation by determining the molecular conden- 



sate fraction after a ramp [3j. More subtly, the multi- 
component condensate is expected to have a complex 
collective mode structure ^3 that will then influence the 
dynamical response through the crossover • The most 
direct signatures are of course in spectroscopy be- 
cause different gap features will correspond to different 
Raman transitions. 

In summary, we have shown that the FR in the ^"^K sys- 
tem involves a three-state Fermi Hamiltonian. Of course, 
while the FM remains only sparsely populated, the char- 
acter of the mean-field ground state shows few qualita- 
tive differences from a single-channel theory, as would 
a Fermi-Bose model in that limit. However, when the 
FM population is significant, the development of weight 
in both (1,3) and (2,3) fractions is revealed in the ap- 
pearance of two length scales in the internal condensate 
wave function. The existence of "Pauli blocking" discrim- 
inates this behavior from that of a Fcrmi-Bosc model. We 
expect signatures of the internal structure of the compos- 
ite wave function will appear in both the collective mode 
response of the condensate and in the dynamics of con- 
densate formation 
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